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We  have  measured  the  ab-plane  thermal  conductivity  (k)  on  single 
crystal  YBa2Cu307-<5  (5  «■  0.08)  in  the  temperature  range  lOK  ^  T  ^ 
300K.  Our  main  focus  is  an  analysis  of  the  slope  change  in  k (T) 
which  occurs  just  below  the  superconducting  transition.  The 
reduced-temperature  Ct=T/TJ  derivative  of  the  normalized  thermal 
conductivity,  d (kVk") /dt  |  t-i,  as  determined  from  the  data,  is 
rather  small,  ^  1.1.  From  measurements  of  the  electrical 
conductivity  on  the  same  specimens  and  application  of  the 
Wiedemann-Franz  Law  we  determine  the  relative  contributions  to 
the  heat  conduction  from  the  carriers  and  the  lattice,  and 
discuss  the  normal-state  phonon  scattering  mechanisms.  Employing 
these  results  we  calculate  the  slope  of  the  lattice  thermal 
conductivity  at  Tc  and  infer  that  the  slope  of  the  carrier 
conductivity  must  be  very  large,  ^  6.  This  result  implies 
strong  coupling  for  some  of  the  carriers. 


PACS  no.'s:  74.70.Vy,  74.30.Ek,  72.15.Eb 


The  thermal  conductivity,  k,  in  the  high-T^  superconductors^ 
is  unique  in  providing  information  on  the  transport  mechanisms 
occurring  in  both  the  normal  and  superconducting  states.  The 
dominance  of  the  phonon  heat  transport  over  that  of  the  carriers 
and  the  importance  of  phonon-carrier  scattering  in  the  cuprates 
are  manifested  as  a  sharp  upturn  in  the  ab-plane  k  for  T  <  T^. 

An  analysis  of  this  feature  offers  the  potential  of  revealing 
information  about  the  superconducting  energy  gap  A  and  the 
strength  of  the  electron-phonon  coupling,  two  fundamental 
parameters  which  are  currently  the  focus  of  much  experimental  and 
theoretical  research. 

In  this  paper  we  analyze  the  slope  change  in  the  ab-plane 
thermal  conductivity  of  two  YBa2Cu307.i  (5  ■>  0.08)  single  crystals. 
Measurements  of  both  the  thermal  and  electrical  conductivities  on 
the  same  specimens  enable  us  to  estimate  the  carrier  and  lattice 
contributions  to  the  total  k.  The  importance  of  phonon-phonon, 
phonon-carrier,  and  phonon-defect  scattering  are  demonstrated. 
Employing  these  results  we  calculate  the  slope-change  in  the 
lattice  conductivity  at  the  superconducting  transition.  We  infer 
that  the  carrier  thermal  conductivity  decreases  sharply  below  T^, 
consistent  with  strong  coupling  for  some  of  the  carriers. 

The  YBajCUsOy.,  (YBCO)  crystals  were  grown  in  zirconia 
crucibles  by  a  self-decanting  flux  method  described  elsewhere,^ 
and  were  subsequently  annealed  in  oxygen  at  450-500 °C.  Both 
specimens  were  high-quality,  mm-sized  (80jim-thick)  single 
crystals  with  p(295K)-250  pQ  cm,  transition  widths  AT^  ^  0.5K, 


and  R=0  at  92K. 
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The  oxygen  deficiency  of  specimen  1  was  estimated  as  &  - 
0.08  from  measurements  of  thermoelectric  power,  performed 
simultaneously  with  the  thermal  conductivity  and  reported 
previously.^  This  value  is  consistent  with  the  value  inferred 
from  the  anomalous  features  in  magnetization  hysteresis  loops, 
which  have  recently  been  shown*  to  correlate  with  c-axis  lattice 
parameter  (a  measure  of  oxygen  content) .  The  magnetization  of 
specimen  2  implies  an  oxygen  content  similar  to  that  of  sample  1. 

The  thermal  conductivity  was  measured  using  a  steady-state 
technique,  employing  a  differential  chromel-constantan 
thermocouple  and  small  resistive  heater  glued  to  the  specimen 
with  varnish.  The  temperature  difference  during  measurement,  AT, 
was  typically  0.3K-1.0K.  Linearity  in  the  AT  response  was 
confirmed  throughout  the  temperature  range  by  varying  the  heater 
power.  Errors  in  k  due  to  heat  losses  are  estimated  to  be  less 
than  2%  except  for  T  >  200K  where  radiation  losses  become  more 
significant.  At  room  temperature  we  estimate  an  error  in  k  of  - 
10%. 

Figure  1  shows  k(T)  for  both  specimens.  Near  room 
temperature  k  -  8-9  W/mK,  which  agrees  with  previous  me-  ^-urements 
on  single-crystal  YBCO.^  The  most  notable  feature  in  chese  data 
is  the  sharp  upturn  and  peak  that  occur  for  T  <  T, .  This 
behavior  is  widely  observed  in  the  cuprates,^  and  arises  from  an 
increase  in  the  phonon  mean  free  path  due  to  the  reduced 
scattering  of  phonons  by  carriers  as  the  latter  condense  into 
superconducting  pairs.  An  analysis  o^  this  slope  change  is  our 
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principal  focus  and  is  discussed  in  detail  below.  We  first 
address  the  behavior  of  k  for  T  >  T,. 

The  contribution  of  the  carriers  to  heat  conduction  may  be 
calculated  v;sing  the  Wi  edemann-Franz  law  (WFL)  ,  which  states  that 
k,/oT  =  L(3,  where  k,  is  the  electronic  component  of  the  thermal 
conductivity,  o  is  the  electrical  conductivity,  and  Lo=2.45xl0‘® 
WQ/K^  is  the  Lorenz  number.  This  yields  an  upper  limit  estimate 
of  K,  and  we  find  at  lOOK,  k^/k  -  0.31  for  both  specimens. 

In  Fig.  2  we  plot  the  lattice  conductivity,  Ki.=k-k,,  where  k, 
is  determined,  as  above,  from  the  measured  electrical  resistivity 
and  the  WFL.  The  slight  upturn  in  Kl  at  T  i  lOOK  for  specimen  2 
is  probably  associated  with  superconducting  fluctuations,  an 
issue  we  address  in  more  detail  elsewhere.®  In  general,  we 
anticipate  three  principal  sources  of  scattering  for  phonons  in  a 
metal:  defects,  free  carriers,  and  other  phonons.  Phonon-carrier 
and  phonon-defect  scattering  yield  a  temperature-independent  or 
weakly  decreasing  with  decreasing  temperature  at  high 
temperatures.^  The  observed  decrease  in  Kl  with  increasing  T  is 
characteristic  of  phonon-phonon  scattering.  The  overall  change 
in  Kl  from  300K  down  to  T^  is,  for  all  specimens,  relatively 
small  (i  20%)  in  comparison  to  that  expected  for  a  defect-free 
insulator  a  factor  of  3).  Thus  with  regard  to  lattice 
conduction,  YBCO  appears  to  be  in  the  complicated  regime  where 
all  three  scattering  mechanisms  are  significant.  These 
observations  are  consistent  with  a  more  quantitative  analysis 
discussed  below. 
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We  now  address  the  change  in  heat  transport  at  the 
superconducting  transition  and  consider  the  ratio  of  the  thermal 
conductivities  in  the  superconducting  (s)  and  normal  (n)  states, 
k®/k",  a  quantity  which  is  most  readily  compared  with  theory.  We 
may  write  the  temperature  derivative  of  this  ratio  as, 


dc 


k^/k") 


(l-Tl) 


dp  ^  dl HK ” 
dt  dt 


(1) 


where  t=T/T;.  is  the  reduced  temperature,  a»K,VK,"/  P"Kl*/Kl'',  and 
ti=k,Vk"  is  the  relative  contribution  of  carriers  to  the  total 
heat  conduction.  The  normal-state  thermal  conductivity,  k",  is 
determined,  for  T  <  T^,  from  the  experimental  data  by  smooth 
extrapolation  of  the  normal-state  data  to  several  degrees  below 
Tj.®  In  Fig.  3  we  plot  k*/k”  versus  t  near  t=l.  The  slopes  are 
d (k*/k") /dt  I  =  --0.4  (sample  1)  and  -1.1  (sample  2)  as  indicated 
by  the  solid  lines.  The  third  term  on  the  right-hand-side  of  (i) 
is  determined  directly  from  the  data  at  T  ^  T^.  This  term  is 
quite  small  and  thus  the  slope  is  determined  by  a  competition 
between  the  first  two  terms  which  have  opposite  sign. 

The  BCS  theory  of  thermal  conductivity  for  weak-coupling 
superconductors®'^  predicts  that,  in  the  limiting  cases  of 
predominant  electron-phonon  and  electron-defect  scattering, 
respectively,  a  is  a  universal  function  of  A/kgT  and  for  A=A®‘^®, 
da/dt|t^.i  =  1.6  (el-ph.)  and  0  (el-def.).  For  strong  coupling  and 
when  electron-defect  scatterings^  ”  predominates,  da/dtl^-i  =  0  as 
in  weak  coupling.  The  case  of  strong  coupling  and  dominant 
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electron-phonon  scattering  is,  however,  much  more  complicated  and 
depends  on  the  form  of  the  electron-phonon  spectral  function, 

(o)  F  (o)  .  For  materials  such  as  lead  and  mercury  where  these 
parameters  are  known  from  tunneling  measurements,  the  very  rapid 
decrease  in  k,  observed  experimentally  is  in  reasonable  accord 
with  theory.  The  difference  from  weak  coupling  arises  both  from 
the  more  rapid  opening  of  the  gap  for  strong  coupling  as  well  as 
the  stronger  frequency  dependence  of  the  quasiparticle  relaxation 
rate.**'^^  The  value  of  da/dt|.,i  =  9  for  lead^®  is  the  largest 
known  for  any  material.  A  slope  da/dt|.,,  substantially  greater 
than  1.6  may  thus  be  considered  evidence  for  strong  electron- 
phonon  coupling. 

The  slope  of  the  lattice  conductivity  for  weak  coupling  and 
predominant  phonon-electron  scattering  can  be  calculated  from  the 
BRT  theory.^®  Tewordt  and  Wolkhausen  have  modified  this  theory 
to  include  the  effects  of  other  scatterers^^  and  crystalline 
anisotropy^®  on  the  lattice  conductivity.  In  addition  to  the 
dependence  on  A/kgT,  the  slope  dp/dtl^.^  depends  on  the  model 
phonon  scattering  rate,  in  the  case  of  YBCO  a  rate  which  must,  as 
discussed  above,  incorporate  the  combined  effects  of  scattering 
by  carriers,  phonons,  and  defects.  For  strong  coupling,  the 
lattice  conductivity  very  close  to  T^  can  be  calculated  from  the 
weak  coupling  theory,  provided  the  BCS  gap  is  scaled  to  the 
appropriate  strong  coupling  value.  ^®  Thus  the  temperature 
dependence  of  P  in  this  regime  is  governed  by  A(T)  which  should 
follow  the  mean-field  expression,  A(T)  -  aCl-t)^'^,  independent  of 
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the  coupling  strength  (a  is  a  scaling  constant  -  1.74A(0)  for  the 
BCS  gap) 

It  is  important  to  note  that  the  change  in  the  lattice 
thermal  conductivity  at  generally  reflects  the  coupling  of 
carriers  to  longitudinal  acoustic  phonons  (presumed  responsible 
for  most  of  the  lattice  heat  conduction)  and  these  phonons  may 
not  be  the  most  important  for  superconductivity.  Recent 
calculations  by  Peacor  et  al.^®  suggest  that  the  coupling  for 
these  modes  in  YBCO  is  rather  weak.  It  is  the  value  of  do/dt|..i 
which  holds  the  most  direct  information  regarding  the  strength  of 
the  carrier-phonon  interaction  relevant  to  superconductivity. 

We  now  calculate  the  lattice  conductivity  slope  from  the  BRT 
theory"^' and  place  limits  on  the  slope  of  the  electronic 
component  for  YBCO  using  our  experimental  results  and  Eg.  (1). 

The  lattice  conductivity  is  written  as,^° 


^8  (  x«e^ 

2k^  v\  j  Jo  (e  *-i)  2 


Jo  2 


(2) 


where  v  =  5,000  m/s  is  the  sound  velocity^^,  0  =  400K  is  the 
Debye  temperature^^,  x=>i(i)/kBT  is  th  reduced  phonon  frequency, 
and  T(x,CjT)  is  the  phonon  relaxation  time.  The  scattering  of 
phonons  by  crystal  boundaries,  point  defects,  carriers,  and  other 
phonons  is  represented  by  ■c‘^(x,T)=  v/d  +  Apd(xT)*  +  A^.xTg  (x,  T)  ( 1- 
+  AppX^T*,  where  d  is  the  crystallite  thickness.  Here, 
g(x,T)  is  the  ratio  of  the  phonon-carrier  scattering  rates  in  the 
normal  and  superconducting  states  as  defined  by  BRT,  equal  to 
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unity  for  T  >  and  a  universal  function  of  A/KbT  for  T  <  T,.. 

The  integration  over  the  variable  C  describes  the  anisotropy  of 
the  electron-phonon  interaction;  it  is  assumed  that  only  the  in¬ 
plane  component  of  phonon  momentum  is  emitted  or  absorbed  by  the 
carriers.^®  Following  reference  17  we  take  where  x  is  a 

constant.  The  coefficients  Ap^,  Ap,,  App  are  determined  self 
consistently  by  numerically  integrating  Eg.  (2)  to  fit  the 
normal-state  data  and  the  data  for  T  ^  20K  where  k  -  k^.  The 
normal-state  curves  for  the  two  specimens  are  shown  as  solid 
lines  in  Fig.  2.  The  slope  dp/dtl^-i  is  then  determined  from  the 
calculated  P(t)  curves.  In  the  inset  of  Fig.  3  we  plot  -dP/dtlp.^ 
versus  the  gap  size  2A(0)/kBTp  for  the  two  parameter  sets.  A 
number  of  experiments  suggest  the  existence  of  two  gaps  in  YBCO, 
associated  with  the  CuOj  planes  and  Cu-0  chains,  with  2Api(0)/kgTc 
-  5-6  and  2  Apj,  ( 0) /ksT,.  -  2. 5-3.5.  The  two-band  expression^®  for  Kl 
is  more  complicated  than  Eq.  (2) ,  but  we  may  place  a  conservative 
lower  limit  on  dp/dtl^^i  by  using  the  smallest  gap  size,  2A(0)/kBT^ 
=  2.5.  This  implies  dP/dt|p.i  -  -4.5  which,  from  Eq.  (1),  yields 
da/dt|t.i  «  6.  A  more  reasonable  picture^®  that  assigns  strong 
coupling  to  the  planes  and  weak  coupling  to  the  chains,  would 
imply  a  plane  contribution,  daPVdtjp.i  -  20.  Our  main  conclusion, 
do/dt|t,.i  ^  6,  indicates  that  strong  coupling  characterizes  some 
of  the  carriers  in  YBCO. 

Since  our  analysis  relies  on  the  calculated  value  of 
dP/dt|t.i,  and  specifically  the  choice  of  scattering  parameters, 
we  have  examined  the  influence  of  the  latter  carefully.  The 
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coefficients  employed  for  sample  1  (2)  imply  that  at  T=100K 
defect  and  carrier  scattering  account,  respectively,  for 
approximately  75%  (60%)  and  18%  (30%)  of  the  lattice  thermal 
resistivity.  Hence  dP/dt|^.i  is  determined  principally  by  the 
relative  weights  of  the  carrier  and  defect  scattering  terms.  We 
find  that  dP/dtl^-i  varies  by  less  than  20%  when  these  terms  are 
varied  within  acceptable  ranges.  We  also  find  similar  results 
when  we  allow  for  a  substantial  deviation  from  the  WFL,  as  would 
be  expected  for  T  <  9  in  a  pure  (electronically  speaking) 
material^  or  a  material  with  very  strong  electron-phonon 
coupling.  Thus  we  believe  that  our  overall  conclusion  regarding 
da/dtlt-i  will  survive  a  more  detailed  and  specific  model. 

As  mentioned  above,  in  the  context  of  the  existing  theory 
values  of  da/dt ( larger  than  1.6  imply  strong  electron-phonon 
coupling  and  predominant  electron-phonon  scattering.  Electron- 
phonon  scattering  predominates  at  only  in  the  purest  of 
conventional  materials  and  one  may  question  whether  it  is 
feasible  for  this  condition  to  hold  in  YBCO.  The  most  obvious 
reason  why  it  may  be  reasonable  is  the  higher  transition 
temperature  of  YBCO;  the  relative  weight  of  phonons  in  the  total 
scattering  of  carriers  is  substantially  enhanced  at  higher 
temperatures.  Secondly,  the  strength  of  the  electron-phonon 
coupling  may  be  substantially  larger  in  YBCO  than  in  any 
conventional  material. 

The  slope  da/dt ( for  strong  coupling  and  predominant 
electron-phonon  scattering^*’  is  directly  proportional  to  the 
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product  of  the  derivative,  d ( A/kg) V^t and  a  frequency  average 
of  the  quasiparticle  relaxation  rate,  r((j).  Employing  the  mean- 
field  expression  for  the  gap  at  T^,  a  gap  size  of  2A(0)/kBT^  -  6 
would  enhance  d  (A/ka)  V^t  |  by  roughly  a  factor  (6/4.1)^  -  2  over 
that  of  lead,  for  example.  It  is  generally  true  that  the  more 
strongly  r((<))  decreases  with  frequency,  the  larger  is  the  slope 
da/dtlt-i*  A  strong  frequency  dependence  can  arise  from  shatj.  ^y 
peaked  structure  in  the  coupling  function,  (o) F (o) .  In  this 
regard  we  note  that  there  is  substantial  evidence  from  Raman 
scattering^'*  that  some  carriers  couple  strongly  to  optical 
phonons,  especially  those  modes  yielding  peaks  in  the  phonon 
density  of  states  F(o)  near  energies  of  15meV  and  45meV.  Recent 
tunneling  results^^  on  Bi2Sr2CaCu20  reveal  such  structure  in  the 
derived  coupling  spectrum  (with  X  -  3.5),  and  similar  results  for 
YBCO  are  anticipated.^®  Thus  it  is  possible  that  uhe 
contribution  to  da/dt|t,.i  from  r{(i))  is  also  large.  As  more 
information  about  0^(0) F(w)  becomes  available  it  will  be  possible 
to  examine  these  hypotheses  in  more  quantitative  fashion. 
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Figure  Captions 

Fig.  1  Thermal  conductivity  versus  temperature  for  two  YBCO 
single  crystals. 

Fig.  2  Lattice  conductivity  (Kl)  versus  temperature,  was 
calculated  by  subtracting,  from  the  data  in  Fig.  1, 
as  determined  from  the  measured  electrical  resistivity 
and  the  WFL.  The  solid  lines  are  calculated  from  Eq, 
(2)  with  A=0  and  parameters  for  sample  1  (2): 
Apd=1.73xl03  (2.80x10^)  Ap,=6.0xl0^  (1.4x10^)  s'^K'\ 

App=40  (19)  s'^K'*  (see  text  for  details)  . 

Fig.  3  The  normalized  thermal  conductivity  near  T^  plotted 

versus  reduced  temperature.  Inset;  values  of  the  the 
reduced-temperature  derivative  of  the  normalized 
lattice  thermal  conductivity,  dP/dt|t,,i,  plotted  versus 
gap  size,  as  calculated  from  Eq.  (2)  using  the  same 
parameters  as  for  the  normal-state  curves  in  Fig.  2. 


